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Abstract 


This  report  contains  the  software  package  for  an 
adaptive,  optimal  suspension  control  system  relative  to 
terrain  random  vibration  disturbances.  The  proposed 
problem  solution  is  shoiim  to  fall  into  two  separate  program 
categories:  a)  recognition  of  the  terrain  and  parameter 
selection,  by  means  of  an  on-board  minicomputer  or  micro- 
processor, b)  optimization  of  suspension  parameters  for 
arbitrary  terrain  configurations  obtained  from  terrain 
statistics  and  executed  on  a centrally*located  stationary 
computer  facility. 

The  interface  between  the  stationary  computer  facility 
and  the  on-board  microprocessor  is  accomplished  by  means 
of  a data  bank  prepared  at  the  stationary  facility  and 
permanently  stored  in  the  meritory  of  the  on-board  micro- 
processor. The  suspension  parameters  are  set  by  a servo- 
control  unit  on  the  vehicle  which  is  activated  by  the  micro- 
processor. The  servo-control  unit  regulates  the  supply  and 
release  of  air  in  the  hydro-pneumatic  suspension  system, 
thereby  increasing  or  decreasing  the  spring  rate  according 
to  the  optimal  requirements.  In  a similar  manner  the 
damper  orifice  size  is  increased  or  diminished  depending 
on  the  required  effective  damping  parameter.  If  need 
arises,  the  vehicle  can  operate  at  fixed  suspension  para- 
meters. The  results  of  the  investigation  are  shown  in 
graph  form. 
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NOTATION 


A(m) 

A(m^) 

a,b 

c (N-sec) 
m 

F(N) 

H(m) 

h(t) 

h, h* 

i, n,m 

k{N) 

m 

K(N) 

m 

«.(m) 
in(kg) 
p (N/m) 
s (m) 
t (sec) 
u (in/sec) 
V(TO^) 

Xf  (in),x^,x^ 
X 

y 


Wave  mumber 

P.S.D.  amplitude  of  road  profile 

Effective  area  of  air  apring 

Subscripts,  axle,  body 

Damping  constant 

Constants 

Force 

Axle  clearance 
Input  function  of  time 

Complex  and  conjugate  complex  funcion  of  time 

Indices 

Spring  rate 

Tire  spring  rate 

Wave  length 

Mass 

Pressure 

Displacement  across  spring 
Time 

Forward  velocity 

Volume  of  air  spring  system 

Displacement  in  direction  of  travel,  also 
vibration  displacements 

Displacement  amplitude 

Relative  vibration  displacement  amplitude 


vi 


Llat  of  Figures 


Figure  la 

Actual  Terrain 

26- 

■27 

lb 

Reconstructed  Terrain 

28-29 

Ic 

Power  Spectral  Density  of  Terrain 

30 

2 

(Xtt)^/g)  vs  u (H  - 0.15m) 

31 

3 

(Xui^/g)  vs  u (H  ■■  0.21m) 

32 

4 

(Xw^/g)  vs  u (H  ■ 0.24m) 

33 

5 

(Xu)^/g)  vs  u (H  ■ o.21m) 

(Constant  Suspension  Parameters) 

34 

6 

Y vs  u (H  ■ 0.21m) 

(Constant  Suspension  Parameter) 

35 

7 

Comparison  of  (Xu  /g)  vs  u (H  - 0.21m) 
between  constant  and  variable  suspension 
parameters 

36 

8 

Block  diagram  of  Suspension  Control 

37 

9a 

K vs  u (H  ■ 0.15m) 

38 

9b 

K vs  u (H  - 0.21ro) 

39 

9c 

K vs  u (H  - 0.24m) 

40 

10a 

C vs  u (H  > 0.15m) 

41 

10b 

C vs  u (H  ■ 0.21m) 

42 

10c 

C vs  u (H  > 0.24m) 

43 

V 


Table  of  Contents 

Abstract  lii 

: List  of  Figures  v 

' Notations  vi 

I 

1 

I 1-  Introduction  1 

1 

2-  Summary  of  Results  4 

3-  Terrain  Analysis  5 

4-  Reconstruction  of  Terrain  Profile  on  Vehicle  7 

5-  Process  of  Averagins  h(t).  Power  Spectral  Density  9 

6-  Design  Constraints  11 

7-  Methodology  13 

8-  Discussion  of  Results  of  the  Investigation  18 

9-  Conclusions  23 

References  24 

Figures  la-lOb  26 

Appendix  A 44 

Appendix  B 60 

Appendix  C 74 

Distribution  List  94 

Report  Documentation  Page  98 


i 


a 

> 

6 

A 

C 

M 

Vi 

P 

T (sec) 
♦ 

uj.RAd 

^Sec 

flRad 

Sec 


•robabllltv  «.ctor  for  4i.pl.ce«nt 

ind  constant  volume 
symbol  for  variation. 

symbol  for  determinant 

Damping  parameter 

Frequency  ratio 

K/k  ratio  of  spring  constants 

m^/rojj 

Weight  factor 
Vibration  period 

PO«.r  .P.ctr.1  4.b.lty  function 
Circular  frequency  Rad/Sec 

Wave  number » (also  a) 


i 


1)  Introduction 

The  investigation  discussed  herein  began  in  1975  when 
the  U.S.  Army  Tank  Automotive  Command,  predecessor  to  TARADCOM, 

I entered  into  a contract  with  Wayne  State  University  to  pursue  a 

study  entitled  "Development  of  Computerized  Vehicle  Suspensions". 

I 

The  investigation  was  completed  in  March  1976.  A final  report  of 

i 

the  work  was  then  submitted  and  approved,  (11.  Subsequently  a new 

f proposal  was  submitted  for  the  purpose  of  establishing  a functional 

algorithm  which  would  allow  the  automatic  control  of  suspension 
I parauneters  on  vehicles  operating  in  rough  terrain.  The  vehicles 

' considered  are  heavy  and  of  the  military  variety;  however,  there  is 

no  inherent  reason  why  the  principles  of  adaptive  control  discussed 
herein  should  not  be  equally  applicable  to  all  vehicles  required 
to  be  operable  at  speed  over  off -road  terrain. 

The  notion  of  automatic  suspension  appeared  first  in  connection 
with  a complete  design  in  a paper  by  Federspiel-Labrosse , [2]. 

I A similar  study  was  carried  out  by  Westinghouse  in  1965,  [3].  The 

I necessity  of  improving  syspension  characteristics  was  soon  put  on  a 

! more  scientific  basis.  Works  by  Bender,  Paul,  Fenoglio,  Karnopp 

^ (at  Massachusetts  Institute  of  Technology  at  the  time  [4, 5, 6, 7]} 

r 

considered  feed-back  automatic  control  systems  for  vehicles.  It 
was  next  shown  that  a transfer  function  can  be  synthesized  (n.  Wiener 
[8]).  Also  Thompson  [9]  considered  the  optimal  active  suspension  for 
• bounce,  pitch  and  roll  control  of  passenger  cars.  None  of  the  above 

! publications  suggested  the  use  of  an  adaptive  control  employing  on- 

i ' board  minicomputer  circuits.  The  reason  for  this  is  that  computer 

i 

j technology  was  not  as  advanced  at  that  time,  as  it  is  today.  Neither 

size  nor  price  form  obstacles  for  the  adoption  of  such  control  methods 
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now,  and  will  be  much  lass  so  in  the  future. 


The  terms  "optimal"  control  and  "adaptive"  control  in 
connection  with  active  vehicle  suspension  designs  require 
further  explanation.  In  a strictly  literary  sense  there  does 
not  exist  a "best"  or  "optimal"  set  of  suspension  parameters 
that  is  realisable,  because  the  "best"  type  of  suspension  is 
one  which  completely  isolates  the  sprung  mass  from  any  form  of 
road  shock.  Any  and  all  suspensions,  no  matter  how  soft,  are 
force  transducers  and  the  softer  the  suspension  the  greater  must 
be  the  excursions  of  the  unsprung  masses,  which  are  not  limitless. 
But,  obviously,  within  given  design  criteria  there  are  better 
performing  suspension  systems  and  the  object  of  this  investigation 
is  to  find  means  to  select  the  best  possible  set  of  suspension 
parameters  that  meet  the  design  criteria. 


The  term  adaptive  control  is  commonly  used  by  control 
engineers  to  indicate  that  the  control  process  is  adapted  to 
the  source  of  the  disturbance  or  perturbation.  Another  alternative 
feedback  control  where  in  the  response  of  the  system  is  compared 
with  a desired  output  and  the  controller  acts  to  minimize  that 
difference  between  the  actual  and  desired  resp>onse. 

In  this  investigation  we  employ  an  adaptive  control  algorithm 
that  allows  pre-sampling  of  the  source  of  disturbance,  namely 
the  terrain  roughness,  from  statistical  data,  (10,11].  From 
the  infinitely  wide  range  of  terrain  configurations  that  exist 
a large  and  substantially  representative  class  of  terrains, 
referred  to  herein  as  model  terrains,  can  be  assembled  and  their 
power  spectral  density  (p.s.d.)  functions  can  be  placed  side-by-side. 
It  can  be  shown  that  in  most  cases  the  distribution  of  the  rough- 
ness follows  similar  patterns  which  can  be  mathematically  expressed 


by  error  functions  (discovered  by  Gauss).  Such  distributions 
are  also  referred  to  as  Gaussian  distributions. 

For  this  investigation  a total  of  13  such  model  terrain 
statistics  were  made  available  and  one  was  selected  for  the 
purpose  of  analysis.  The  conqputer  algorithm  dedicated  to  the 
on-board  microproce'^r^or  was  simulated  and  combined  with  the 
off-board  terrain  and  optimal  suspension  parameter  selection 
process,  so  as  to  close  the  look  of  the  process.  The  following 
chapters  will  explain  in  detail  the  logic  of  the  program  and 
its  operational  features. 
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2)  Summary  of  Results 

a)  Raconstruct ion  of  terrain. 

The  reconstruction  of  the  terrain  by  means  of  a 
fast  Fourier  analysis  and  regeneration  of  measurable  accelera- 
tion values  at  vehicle  unsprung  and  sprung  mass  yielded  very 
accurate  results.  The  terrain  as  described  in  the  form 
of  tabulated  data,  (101  is  plotted  in  Figure  19.  Of  course,  the 
accelerometers  employed  must  have  a nearly  linear  output 
relative  to  the  input  within  the  frequency  domain  of  their  use. 

b)  Central  Stationary  Computer  Simulations 

The  program  for  the  computation  of  the  Optimal  parameter 
matrix  has  been  written,  documented  and  executed.  The  results 

have  been  plotted  and  compared  with  fixed  suspension  parameter 
systems.  At  comparable  vertical  acceleration  levels  the 
realizable  speeds  (speeds  made  good)  of  systems  with  adaptive 
suspension  control  exceeded  those  speeds,  obtainable  with  fixed 
suspension  parameter  systems  (passive  suspension) . If  vehicles 
are  operated  at  the  recommended  speeds  the  vertical  accelerations 
of  the  sprung  masses  are  expected  to  be  sufficiently  reduced 
to  extend  the  life  of  vehicle  components.  But  without  statistically 
significant  evidence  no  quantitative  statement  can  be  made. 
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3)  Terrain  Analysis 


r 


The  signifiance  of  terrain  configurations  is  that  once 
their  pattern  is  estaiblished  it  is  possible  to  reconstruct 
them  without  a substantial  loss  of  accuracy  relative  to  the 
response  characteristics  of  vehicles  passing  over  such  terrains. 

From  terrain  measurements  (approximately  120  for  each  terrain 
sampled  over  a length  of  circa  120  meters)  one  can  obtain 
equivalent  Fourier  series  which  are  considered  as  representa- 
tive for  each  terrain  configuration,  (see  Figure  1) . 

The  computer  program  for  the  "Fast  Fourier  Transform" 

[12]  was  successfully  executed  and  the  coefficients  of  each 
sine  and  cosine  term  were  obtained.  The  reconstruction  of  the 
terrain  from  the  128  terms  considered  in  the  series  showed  an 
accuracy  of  between  97%  and  99%  relative  to  distribution  of 
the  measured  data. 

Since  random  data  do  not  repeat  themselves,  it  is  necessary 
to  take  assembly  averages.  All  terrains  considered  herein  are 
considered  to  be  ergodic,  that  is  to  say,  that  different 
samples  of  the  same  terrain  when  averaged  over  the  length  of 
the  sample  yield  essentially  the  same  value.  Taking  the 
squared  values  of  the  measured  elevations  from  a mean  level 
(squaring  is  necessary  to  avoid  getting  simply  the  algebraic 
mean  value)  and  dividing  the  sum  of  all  squared  elevations 
by  the  number  of  measurements  one  obtains  the  mean  square 
elevation  of  the  terrain  profile  and  the  root  mean  square 
value,  respectively.  Hence  for  an  ergodic  input  the  averages 

l; 

j 

I 
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of  all  samples  must  be  nearly  alike  for  a specific  type 
of  terrain.  Note  that  we  do  not  assume  that  terrains 
are  stationary,  which  would  imply  that  the  road  elevations 
repeat  themselves  within  multiples  of  a base  length  of  the 
terrain,  because  terrains  do  change  over  longer  distances. 

In  fact,  it  would  be  necessary  to  collect  an  infinitely  large 
number  of  samples  to  characterize  all  types  of  terrain. 

Prom  a practical  point  of  view,  however,  it  is  possible  to 
select  a finite  set  of  terrain  configurations  which  a parti- 
cular vehicle  will  encounter  during  operation  and  to  determine 
the  most  suitable  form  of  suspension  characteristics  for 
operation  in  such  terrains.  The  set  of  terrains  so  selected 
are  referred  to  as  model  terrains. 

It  is  assumed  that  the  wheels  of  the  vehicle  passing 
over  the  terrain  will  stay  in  contact  with  the  terrain  profile 
at  all  times.  Also,  since  the  total  length  of  the  measured 
sample  extends  over  100  times  the  average  distance  between 
value  sets  the  mean  square  value  of  elevation  for  each  wheel 
may  be  assumed  to  be  the  same.  In  contrast  to  this  assumption 
if  one  were  to  analyze  a specific  wave  disturbance  of  a 
fixed  wave  length  the  wheel  elevation  of  each  wheel  would  be 
different  at  any  one  time  and  It  therefore  would  become 
necessary  to  consider  the  phase  relationship  between  each 
wheel  (the  distinction  is  made  to  indicate  the  difference 
between  a stochastic  and  a deterministic  input )» 


4)  Reconstruction  of  Terrain  Profile  on  Vehicle 


The  vehicle  is  to  be  equipped  with  low  frequency 
sensitive  accelerometers  on  the  superstructure  (body)  and 
high  frequency  accelerometers  on  the  axle  (unsprung  mass). 
With  lOjj  the  pro-rated  sprung  mass,  the  pro-rated  unsprung 
mass,  )c  the  suspension  stiffness  rate,  K the  tire  stiffness 
rate  and  c the  shoc)(  absorber  constant  we  can  write  the 
equations  of  motion  in  the  form 

«»b  ^b  *b  ^ ■ •'’‘a  ■ °*a  “ ° 

(1) 

m-il-  + (K+)c)x_  + ox.  -kx.  -cx.-  Kh(t) 
a a a a D o 

The  displacements  in  (1)  are  x.  (body)  and  x (axle). 

D a 

h(t)  is  the  wheel  lift  and  drop  due  to  terrain  roughness. 
Adding  the  two  equations  (1)  we  obtain 

Vb  ^ *a*a  ^ ^'’‘a  " (2) 

Xa  “ if 


From  (2)  we  obtain  the  terrain  profile  h(t),  namely 


h(t) 


3c. 


(\) 


+ 3t,("’a)  + X, 

a jr  a 


(2a) 


If  X.  and  X.  are  measured  and  converted  into  electronic 
D a 

impulses  then  (3)  can  be  obtained  by  integrating  twice  and 
inverting  the  sign  on  sumning  amplifiers.  Alternately,  the 
accelerometer  outputs  could  be  digitised  and  the  integration  done 
numerically  on  the  microprocessor.  In  either  case  the  the  result  of  (2a) 


yields  a reading  of  the  reconstructed  terrain  profile. 

In  order  to  demonstrate  the  accuracy  of  the  reconstruction 
process  the  investigators  converted  the  spatial  Fourier  series 
of  the  terrain  into  a time  series  based  on  a constant  speed 
motion  of  the  vehicle  as  follows. 

Let  the  n^*'  term  of  the  series  be  sin  (nax^)  where  n is  an 
integer,  a is  the  wave  parameter  (a"2n/t)  where  « is  wave  length 
in  meten)  and  is  the  forward  motion  displacement.  Then 

x^-u.t  (4) 

where  u is  speed  and  t la  time  the  argument 

naxj  »n(2Ti  u/t)‘t  ■ n-wt  (5) 

where  in  (5)  u)  ■ 2nu/t  is  the  circular  frequency  of  the 

harmonic  of  ”?*'  meter  wave  length.  Referring  back  to  (1) 

we  can  compute  x^,  x^  relative  to  each  and  every  term 

contained  in  the  Fourier  series  representing  h(t).  These 

are,  in  fact,  the  signals  sensed  by  the  on-board  accelerometers. 

X.  and  X is  then  reintroduced  into  (2a)  to  obtain  h(t)  and 
D a 

the  output  is  compared  with  the  input.  The  correlation  proved 
to  be  on  the  order  of  a fraction  of  a percent. 
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5)  Process  of  Averaging  h(t)y  Power  Spectral  Density. 


The  terrain  elevation  may  be  expressed  in  terms  of 


discrete  (measured)  values  or  as  a continous  function  of 

.th 


distance  or  time.  The  n term  of  a Fourier  series  ex- 
pansion of  the  terrain  function  may  have  the  form 


h(t)„  = h^^  cos  nut  + h^„  sin  nut 
n nc  ns 


(6) 


where  h__and  h^„  are  the  coefficients  associated  with  the 
nc  ns 


th 


sine  and  cosine  terms  of  the  n order.  Hence 


2 2 2 2 2 
h (t)^  = h cos  nut  + h „ sin  nut  + h h sin  2nut  (7) 
n nc  ns  nc  ns 


Averaging  the  above  with  respect  to  one  period  of  length  t 
we  obtain 


h^  f 2 

nc  I cos  nut  dt  + 


sin  nutdt 


lls  1 

T 

^nc  ^8  I sin, 2 nutdt 

T 


(8) 


Now  the  last  integral  on  the  left  side  of  (8)  integrated 
over  the  full  period  x vanishes  and  (8)  becomes 


-7  ^ k2 

h^  = ( nc  + ns  ) = ( nc  ns  ) 
2 2 2 


(9) 


Each  of  the  squared  amplitude  values  h belongs  to  a discrete 


frequency  of  order  n . Then, the  average  of  all  such  values 
over  the  frequency  spectrum  considered  herein  is 


6)  Design  Constraints 

The  axle  clearance  is  usually  determined  by  several 
design  criteria  other  than  softness  of  ride.  Therefore,  it 
is  necessary  to  avoid  the  incidence  of  the  axle  striking  the 
axle  stops.  The  probability  that  the  maximal  displacement  of 
the  axle  relative  to  the  body  Y will  lie  between  ± ay  = 4 H 
where  a is  a number,  H is  the  available  axle  clearance  and 
y is  the  expected  suspension  deflection, is 

Prob  [-oy  — y{t)  £ ay]  ■ j ” e dy/y  (2ir) 

J-ay 

(14) 

and  for  a * 1,2,3  this  probability  is  [14] 

a = 1 Prob.  = 68.3% 

a = 2 Prob.  = 95.4%  — Gaussian  or  normal 

distribution 

a ■ 3 Prob.  = 99.7% 

This  means  for  a » 3, in  only  3 out  of  1000  working  cycles 
would  Y exceed  the>  allowable  axle  clearance  H and  therefore 
strike  the  axle  stops. 

In  seeking  an  optimal  set  of  suspension  parameters , 
design  constraints  imposed  on  the  system  limit  the 
average  maximal  displacement  across  the  suspension  y to  1/u  = 
1/3  of  available  axle  clearance  H.  For  each  terrain  profile 
there  exists  a set  of  suspension  parameters  k and  c that 
minimizes  the  maximal  body  accelerations  |X|^|  at  a certain 
average  forward  speed  u and  at  the  same  time  satisfies  the 
design  constraint  that  the  displacement  across  the  suspension 
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should  be  within  the  limit  H 99.7%  of  the  time. 


At 


other  speeds, u, the re  are  sets  of  suspension  parameters 
which  are  not  optimal,  yet  satisfy  the  design  constraint. 

These  parameters  are  optimal  to  the  extent  that  the  maximal 
displacement  across  the  suspension  elements  is  utilized. 

Since,  in  general,  the  acceleration  transmissibility  depends 
inversely  on  the  deflection  across  the  suspension  and  directly 
on  the  forward  speed  of  the  vehicle, it  is  advantageous  to 
exhaust  all  the  available  axle  clearance  (without  frequently 
striking  the  axle  stops) , so  as  to  obtain  the  sofest  ride 
possible.  In  Figures  2,3,4  we  show  the  acceleration  in  units 
of  "g"  of  the  vehicle  body  versus  "speed  made  good"  for  three 
differently  chosen  axle  clearances  H ■ .15m,. 21m,  .24m.  All 
of  these  represent  optimal  data,  as  explained  above,  relative  tc 

variable  suspension  spring  rates  and  shock  absorber  constants. 
Figure  5 portrays  the  acceleration  transmissibility  versus 
"speed  made  good"  for  constant  spring  rate  and  damping  constant 
(passive  suspension) . In  all  cases  the  results  obtained  refer 
only  to  one  terrain  profile  (see  Appendix  1) . The  formulas  used 
to  compute  the  set  of  optimal  spring  rates  and  damping  constants 
are  presented  in  Appendix  2.  For  the  purpose  of  showing  the 
effects  of  speed  on  suspension  deflection  we  show  in  Figure  6 the 
relative  axle  displacement  y vs.  u for  fixed  spring  and  damping 
parameters.  Figure  7 shows  comparison  of  acceleration  trans- 
missibility. The  meaning  of  the  data  is  discussed  on  pane  18. 


Methodology 


The  control  process  is  schematically  shown  in  Figure  8 
(Figure  22  of  (11)  in  the  form  of  a block  diagram.  At  the  right 
hand  lower  corner  are  shown  the  off>board  logistic  operations 
executed  on  any  stationary  computer  installation.  Given  a 
discrete  number  of  different  model  terrain  statistics  and 
the  descriptions  of  the  vehicle  masses  and  tire  or  track 
stiffness  optimal  parameter  matrices  (see  Table  1)  are 
obtained  which  give  the  best  possible  combinations  of  spring 
rates,  damping  constants  and  corresponding  forward  speeds. 

The  data  are  placed  in  memory  storage  on  the  on-board  micro- 
processor. 

The  computation  of  optimal  parameters  requires  many 
programming  steps.  These  are  documented  and  attached  to 
this  report  in  Appendix  2.  Speaking  in  general  terms,  the  mean 
square  terrain  profile  amplitude  is  amplified  or  attenuated 
depending  on  spring  rate  k,  damping  constant  c,  and  forward 
speed  u.  Only  two  displacement  values  are  of  significance. 

These  are  a)  the  displacement  across  the  suspension  as  mentioned 
in  the  foregoing  paragraph  and  b)  the  absolute  displacement  of 
the  sprung  mass  because  the  peak  accelerations  are  directly 
proportional  to  that  displacement. 

Thus : 

peak  acceleration  * frequency  squared  times  displacement  or 

l*cl  -"^*1^1 


13 


and  (14)  written  in  units  of  g (gravitational  constant) 

becomes  |i  | - |x.  | (14a) 

^ q ^ 


It  can  be  shown  (see  appendix  2)  that  |X|^|  is  a function  of 
I Y I neunely 

UkI  • I*! 

V (15) 

where 

^ " c/2'1  ; 

2 2 2 2 2 
u)  - 4it'^u  /«  ; u)^  - )c/mjj  . 


I Instead  of  Y we  use  the  permissible  average  relative 

displacement  amplitude 

y = H/a  « H/3. (n-3) 

The  mean  square  terrain  elevation  amplitude  is: 

R » 1 E h^  ; (m  - 120) . 

1 Im  1 " 


(16) 


(17) 


-2  -2  -2  2 

The  ratio  y /h  , where  y » H /9,  is  a specific 
number  deptending  only  on  the  axle  clearance  H and  the 
measured  terrain  roughness  expressed  by  the  Fourier  series 
amplitude  values  h^.  In  appendix  2 it  is  shown  that: 


-2  -2 
(yVh'^)  « 


(18) 


"•b 

a 


l+l]^+r.^n^ 


r 2A 
(n 


% 

-I  )-l 

ir 

a 


2 

where  in  (18)  U ■ K/m  is  a Itnown  parameter.  Solving  (18) 

A A 

2 

for  n as  a function  of  k and  c we  can  infer  the  forward 
speeds  u^  pertaining  to  the  set  of  solutions,  since 


14 


or 


(19) 


Y)r^/(2itymJj) 

Substitution  of  the  solution  into  (15)  yields  the  absolute 
expected  displacement  amplitude  which  should  be  minimized 
for  best  possible  ride  quality.  Again,  assuming  ^constant 
one  can  find  a unique  set  of  values  for  u and  k satisfying 
this  requirement  and  after  further  iteration  by  varying  C 
(or  damping  constant  c)  the  best  parameter  configuration  u^, 
k^,  c^  emerges  whereby  the  subscript  o denotes  optimum.  It 
follows  that  for  each  terrain  only  one  forward  speed  will 
yield  optimal  performance.  For  other  speeds  there  is  a 
parameter  set  k (u) , c (u) , different  from  the  optimal  values, 
that  renders  (15)  a relative  minimum.  The  above  described 
procedure  of  optimization  has  been  coded  in  Fortran  language 
and,  as  far  as  the  principal  investigator  is  aware,  it  is 
orginial. 

Prior  to  the  development  of  the  above  described  procedure 
the  investigators  attempted  to  employ  the  optimization 
method  of  Bender  et  al  of  the  Massachusetts  Institute  of 
Technology.  According  to  that  method  optimal  system  responses 
are  obt'a'ined  by  minimizing  a so-called  penalty  function  P which 
is  equal  to  the  sum  of  the  absolute  displacement  and  the 
relative  displacement  across  the  suspension.  Since  the  re- 
lative displacements  depend  not  only  on  spring  rate  and 
damping  but  foremost  on  the  available  axle  clearance  one  can 
develop  a so-called  trade-off  curve  of  acceleration  trans- 
missibility  versus  mean  terrain  amplitude  and  available  axle 
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clearance.  Obviously,  the  greater  the  axle  clearance  the 
less  is  the  trasmitted  acceleration,  or  for  a fixed  acceptable 
acceleration  transmissibility , the  greater  is  the  vehicle 
speed  with  which  it  can  move  over  the  rough  terrain.  If 

> - -44  (20) 

P * px^  * <xy  , (p  is  a number,  10  <p<10  ) 

then 

dx  ■“ 

dP  = p b ♦ a dy  =>  o 

air  (21) 

^ = c ^ + a ^ = o 

dc  dc  He 

yield  values  of  k and  c that  minimize,  or  maximize  (20) . 

The  differentiation  (21)  yields  a very  complicated  algebraic 
polynomial  of  the  eight  power  in  the  variables  and  even 
though  the  equations  are  coded  for  data  processing  the  results 
obtained  cannot  be  readily  verified.  They  indicated  that  the 

penalty  function  could  be  represented  by  a relatively  flat 
surface  and  therefore  is  insensitive  to  small  variations  in  the 
coefficients  )c  and  c on  which  Xj^  and  y depend  implicity.  The 
results  yield  one  set  of  optimal  parameters  and  c^,  practically, 
f independent  of  variations  in  the  parameter  c.  Thus,  this  approach 

I 

was  later  discarded. 

To  implement  the  results  of  the  optimization  developed  in  this 
Study,  a Table  1,  prepared  especially  for  the  vehicle,  is  then 
placed  together  with  the  mean  square  terrain  amplitudes  in  the 
on-board  computer.  The  computer  is  also  programed  to  reconstruct 
terrain  profiles  and  to  obtain  the  mean  square  amplitude  value 
of  the  terrain  over  which  the  vehicle  passes.  It  then  selects 
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t 
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the  nearest  amplitude  value  model  terrain  and  the  nearest 
speed  for  which  pre-calculated  parameter  sets  are  available. 

The  parameter  set  (k^  and  c^)  information  is  relayed  to  the 
servo  control  unit  for  implementation.  Then,  the  servo  control 
unit  changes  inflation  pressures  of  the  air  spring  and  orifice 
size  of  the  shock  absorl^erf  to  the  desired  levels. 

In  order  to  execute  the  control  function  a comparison 
between  instantaneous  and  recommended  suspension  parameter 
values  is  made  based  upon  the  measured  forward  sp>eed.  The 
optimal  speed  will  be  displaced  in  sight  of  the  vehicle 
operator.  Then,  the  servo-control  unit  will  change  inflation 
pressure  of  the  air  bellows  springs  and  orifice  of  shock 
absorbers  to  the  desired  level.  Since  the  terrain  profile 
will  be  continually  scanned  within  a predetermined  length 
(say  100  meters) , the  operation  constitutes  a closed  loop 
adaptive  control  process. 
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It  is  apparent  that  the  economic  benefits  of  optimal 


suspension  control  reside  in  either  extended  periods  of 

operation  between  overhauls,  or  in  reduced  time  to  cover 

distances  (higher  average  speed)  or  in  both  of  these 

advantages.  It  is  also  clear  that  performance  improvements  require 

an  initial  investment  (no  matter  how  small  in  relation  to 

the  total  equipment  cost) . Whether  or  not  the  expected 

benefits  outweigh  the  cost  is  a question  of  judgement. 

But  the  gain  in  speed  over  a conventional  passive  suspension 
system  can  be  demonstrated.  Let  terrain  profile,  sprung 
and  unsprung  masses,  tire  spring  rate  and  axle  clearance  be 
equal  for  two  vehicles,  one  with  passive  and  one  with  active 
suspension.  Figure  7 clearly  shows  that  the  acceleration  trans- 
missibility  through  the  passive  suspension  system  is  greater 
than  that  of  the  active  suspension  at  comparable  forward 
speeds. 

Transmissibility  ratios  are  never  unique.  By  this  we 
mean  that  bi-quadratic  equations  can  yield  two  distinct 
real  roots  for  each  circular  frequency  value  or  forward 
speed  considered.  Frequency  spectra  have  maxima  at  a 
finite  frequency  value  and  minimum  at  o and  °°.  On  either 
side  of  the  frequency  giving  the  maximum  lie  higher  and 
lower  frequency  values  yielding  the  same  transmissibility. 

For  reasons  other  than  transmissibility  one  cannot  ma)ce 
springs  arbitrarily  soft  and  dashpots  ineffective.  The 
parameter  search  is  then  restricted  to  the  frequency  domain 

that  yields  reasonable  suspension  parameters  at  practically 
realizable  forward  speeds. 
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In  figures,  9,  10  we  show  the  effects  of  parameter 

variation  k and  c, which  is  the  dimensionless  parameter  of 

the  damping  coefficient  c,on  speed  u for  the  optimized 

suspension.  The  spring  force  variations  may  be  realized 

in  two  ways.  Since  the  rate  is  defined  as  [15] 

k ■ d (pA)  * to  A+p  dA  (22) 

3s  3s  3? 

one  has  the  choice  to  alter  either  the  effective  area  A 
relative  to  stroke  s and  let  the  pressure  be  nearly  constant 
or  to  vary  the  pressure  relative  to  stroke  and  let  the 
effective  area  be  constant.  The  former  of  the  two  choices  is 
realized  by  shaping  the  air  cushion  requiring 

dA  ^ k (m^/m)  (23) 

Si”  p 

The  latter  implies  that 

to  = k (N/m^)  (24) 

3s  A 

However,  the  spring  rates  are  non-monotonia lly  increasing 

or  decreasing  function  of  u giving  the  data  listed  in  Table  I. 

From  the  recorded  forward  speed  and  the  recorded, reconstructed  , 

-2 

terrain  function  yielding  the  mean  square  value  h the  appropriate 
suspension  parameters  are  selected.  The  amended  data  matrix 
(Table  1)  will  yield  body  acceleration  values  which  are  relative- 
ly best  for  the  recorded  forward  speeds. 


Table  I 


H ■ .15m 

u 

c 

8 

1000 

1132 

12 

3000 

980 

16 

6000 

462 

20 

9000 

1132 

24 

14000 

706 

28 

18000 

800 

m/sec 

N/ra 

N aec/m 

H - .21m 


24  17000 


28  24000 


H * .24m 

u 

ds 

c 

8 2000 


12  5000 


16  9000 


20  14000 


24  11000  3129 


28  28000 


N sec  ra/sec  N/m 
m 


At  axle  clearance  H = .15m 


= 16  m/sec. 


k » 3000  (N/m) , c = 1132  N sec 

m 


are  optimal  values.  So  are 

u > 8m/sec, 

c ■ 267  N-sec 
m 

for  H » .21m 

and  u = 12m/sec, 

c ■>  422  N sec 


= 2000  N/m 


k = 5000  N , 
m 


H - .24m  . 


The  improvement  can  be  deduced  from  the  comparison  between 
fixed  (assumed)  suspension  parameters  and  optimal  parameter 
response* Figure 7.  In  closing  it  shall  be  mentioned  that 
remarkable  advantages  can  be  realized  with  the  proposed 
suspension  control  system. 

1)  Recommended  speeds  for  optimal  control  can  be  posted 
at  the  dashboard  (digital  read-out) 

2)  For  whatever  reason  the  operator  deems  necessary  he 
c-nn  switch  to  passive  control  or  manual  control,  employing 
the  control  scheme  based  on  (24) , Table  I. 

Example 

A vehicle  negotiates  a terrain  described  in  [10]  and 
displayed  in  Figure  la.  The  terrain  is  Fourier  transformed 
and  axle  and  body  accelerations  are  monitored  for  a vehicle 
forward  speed  of  21m/sec  based  on  the  assumed  data  given  in 
Table  2.  The  terrain  roughness  is  now  reconstructed  by  direct 
integration,  namely, 

h(t)  = ma  X + (x.dt)dt+x.  ”*b  (2a) 

-JT  * * -JT 

giving  terrain  figure  lb.  The  mean  square  value  (10),  (17) 

are  then  introduced  into  (18)  where  the  permissible  mean 

_2 

square  suspension  deflection  amplitude  y is  given  1/9  of  the 

2 

available  square  axle  clearance  H e.g. 

H » .21m,  -=  .0442m^  H^/0  = .00049  = y^  ; 

_ 2 
y = .07  in.  For  each  set  of  parameters  k,  n that 


satisfies  (18)  there  exists  a forward  speed  U)  that  satisfies 
(19)  . 


T^ble  2 

Technical  Data 


"■  ' Tandem  Axle 

Sprung  Mass  M:  4.087  x 10^  )tg  per  wheel  set 

Unsprung  Mass  m:  3.04  3 x 10  )tg  per  wheel  set 

Sprung  Mass 

2 2 

Natural  frequency  60  cpm  » 1 Hz;  * 39.48/sec 

2 

Spring  rate  * 161,354  kg/sec 

2 

Spring  rate  K » 1,201,376  kg/sec 

Damping  Constant  c = 17,976  kg/sec 

1/2 

Damping  Parameter  f.  * c/2  (kM)  « .7 
Wave  Length  L * 100m  (m  ■ meter) 

Speed  u * 25  m/sec  (maximum),  20,  15,  10,  5 m/sec 

Double  amplitude  of  fundamental  wave  of  length  100m:  .1  m (10cm) 

Axle  clearance:  (Static-Source) 
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Any  person  who  is  called  upon  to  review  a prospective 
design  for  possible  adoption  in  the  course  of  product  planning, 
must  ask  whether  the  benefits  gained  outweigh  the  costs  of  the 
proposed  suspension  control.  This  judgement  will  depend  on  many 
factors  most  of  which  are  unknown  to  the  author  of  this  report  at 
this  point  in  time. 

In  connection  with  such  considerations  it  is  fair  to  assess 

the  expected  performance  of  the  system.  Current  expectations  are 

that  best  results  will  be  obtained  in  terrains  for  which  power 

2 2 

spectral  desitles  vary  between  . Im  /rad/m  to  .0001m  /rad/m  in  the 
frequency  range  from  2 to  6 rad/meter.  In  terrains  of  substantially 
greater  roughness  manual  control  may  be  preferable  and  in  terrains 
which  are  much  smoother  fixed  suspension  control  may  be  more 
advantageous.  Both  modes  could  be  made  available  in  this  control 
system.  Also  note  that  no  changes  in  the  parameter  settings  of  the 
suspension  will  occur  unless  speed  or  terrain  configurations  or  both 
change  more  than  set  threshold  levels  during  the  operation  of  the 
vehicle. 

Actual  perfromance  Increases  to  be  expected  would  have  to  be 
determined  at  the  time  that  design  trade  offs  are  being  made. 
However,  based  on  the  simplicity  of  the  system  as  compared  to 
others  being  contemplated,  an  adaptive  suspension  control  mechanism 
appears  to  offer  the  most  attractive  means  of  optimizing  ride 
control  in  military  vehicles. 
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Actual  Terrain  Profile  Figure  1(a) 
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POUER  SPECTIML  DENSITY.  n<n/(  RAOIRNS/m 


Powar  Spactral  Danslty  of  Torrain  Figura  1(c) 
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BLOCK  diagram  OF  ADAPTIVE  CONTROL 
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appendix  a 


Appendix  A 


A program  has  been  established  containing  the  program 
steps  amd  documentation  as  attached. 

The  basic  feature  of  the  program  is  the  subroutine  "Fast 
Fourier  Transform"  (Reference  the  Fast  Fourier  Transform  by 
E.O.  Brigham,  Prentice  Hall,  1974,  pp.  148-171).  The  program 
has  been  transcribed  using  complex  variables.  It  serves  the 
purpose  of  resolution  of  the  recorded  input  signals  into  its 
harmonic  components  within  a time  interval  -T<t<T  wherein  T 
is  the  half  period  as  considered  being  sufficient  in  duration 
to  characterize  the  random  signal  as  that  pertaining  to  a 
stationary  process.  Next,  with  the  help  of  the  system 
equations  (1)  the  perturbation  signal  h(t)  is  regenerated  and 
the  PSD-curve  developed.  This  newly  generated  ensemble  is 
then  compared  with  one  that  is  based  on  the  actual  terrain 
configuration. 
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VESUS 

START 


) 


I 


\ 

COMMON  VEDATA,  PROF I L 


CALL  ROUGH (MP,SH) 


CALL  PSDCAL(NP,SH,S1) 


CALL  FFT(SH,NP) 


' CALL  XANXB (SH,XA,DX2A,DX2B,NP) 


// 


-J 

HCAL(DX2^,DX2B,XA,SH,NEW,NP^•■ 


//  CALL  PSDCAL(NP,SHNEW,S2 


-t 


COMPARE  SI  & S2 


WRITE  S1,S2,  DIFF  IN  SI  & S2 


(RETURN  N 
STOP  J 
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SU9RCUTINE  FrTIA»N»NRI 

O'^O’ 

CC^RIEX  AINRI  »U»H»T»C'*<>LX 

C 

OIVICING  AU  ELE'^ENTS  BY  NR 

C^O? 

OC  \ J«1*NB 

t 

AIJ)  • AIJI/NB 

C 

REORDERING  THE  SEQUENCE 

-rjos 

NB02  - NBFZ 

CO  )6 

NBH  • NS  - 1 

«''07 

J - 1 

oo* 

OC  * 1 ■ l»NBNl 

mo 

IF(L.GE.J)GD  TO  2 

'OK 

T ■ AU) 

0011 

AU)  • AID 

''O'  7 

A<1  ) « T 

roi  0 

> 

K«NRD2 

c 

WRITE(6*lO))(A(I)*l«l»6) 

''“.U 

lol 

FfRNATUFlS.B) 

'•oi^ 

y 

rFIK.GEU)G3  TO  4 

J - ) - K 

•'1’ 

< • KF? 

GO  T : 3 

1 

h 

J - J^K 

c 

CO^ALTAT icn  CF  fft 

• • ;• 

Cl  • :.145??65R59Q?93 

C-  - •>  t 

DC  6 *'■1  *N 

- > ^ 

1'  • ( I.'' #0.0  1 

" > ■*  " 

wC  a >aa V 

K a *»E/? 

w • r'*Pix(f':s<«*i/K)»  -siN)si/x)) 

C'T  b Jal.K 

- 1 » 

:r  5 L» 

L =•<  « L^K 

• -• 

* « ^UPK  ) 'J 

-'>  1 ' 

-)LPK)  • A(L)  - T 

r 

wt:  !TC(  rtlC2)  AUP*') 

• ' . 

K’ 

?F1  ■'.O  ) 

1 ' 

i<L  ) « A<L 

"J-.  • 

*1 

L'  » 1.-*  W 

■JrTJ^N 

■ ' 1 A 

.•  r=F£f’»  IO.P'^C:IC»SjU^lCr»^CL  IS’.K'Oc:'<»L..i:.‘.r-l' 

• • r c c 

jf-a  M4MC  a per  , tlNfCN’  • R? 

• • . 

SCO-iCc  STaTFN=N''S  • 35.PP?'‘.A'A*'  Sl'f  • 

• 1 

• NC  r f : vjcs'ics  36Nc»«at=: 


■'Efi'^lNAL  system  F0'»T«!AN  CI41336) 


main 


-'301 

OTQT 

0*115 
'*'»06 
r •).;-» 

cinp 

r‘)''Q 

• C*''  I 0M5 

• S’-" I S’ 
•statist 


c 

c MCAL  starts  hER  £»*••••  *•**•> 

C IF  tlE  400  THE  Two  EOU4TICN3  QF  “OTION  FCP  B30Y  AND  AXIE 

C THEN  WE  GET  REQUIRED  RELATION  TO  CB’AlN  TERRAIN  ROUGHNESS 

SUBROUTINE  HCAL(0X2A*0X2B>XA.H»NP| 

COMPLEX  H(NO|«xA(NP)#DX2m(N»|,CX2R(NP) 

C0PM:NTVe0AT4/AH»BM.TK.SY*C.PI 

Cl  « AM/TX 

C?  • RM/TK 

DO  1 I • I*NP 

1 HID  • XAII)  ♦ Cl*0X2A(I»  ♦ C?*0X2«(n 
RETURN 
ENO 

IN  EFFECT*  IO»ERCDlC*SCURCE.NOL  IS’»NOOErx»LCAO. nomap 
IN  EFFECT*  NAME  • HCAL  • LINECN’  • 57 


ICS*  SuURC:  statements  • 

ICS-  NC  diagnostics  generated 


9,PR03FA*  SIZ:  - 
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'llf.'l  I".  A^ 

*t  c**!  NAL 

srs"''?'*  cqrtRAN  G<a1336>  IAIN  r’-'>6-T 

C 

C 

FUNCTION  SU3PR0CRAN  FQR  aTH  OROFR  RJNSA-RUTTA  '•F'^HOD 

PUNCTirN  RUNGEIN»T»F,X,HJ 

INTEGER  RUNGE 

01  »»EFSinN  PHI  <5P1  .SAVEYISCUTIN)  .FIN) 

^Tn«, 

OA'A  H/0/ 

C 

»«•••♦! 

r:.:<s 

GC  TC  11. 2.3.  A. 5).** 

'“O'.? 

1 

RUNGk«l 

^ i 

RETUFN 

2 

DC  22  J-l.N 

ic 

SAVEYIJ)  • YIJ) 

'•■>11 

BHI  IJI  • FU) 

'•'11  ■» 

22 

rCJ)  • SAVEYIJI  ♦ .5*H*FIJ) 

'■'>17 

A « y A .5^H 

C31  A 

RUNG!  • 1 

r.  • 1 ^ 

RE''Ul>N 

''7U 

3’  J-l.K 

'■71' 

PHU  J)  • PHI  ( J)  ♦ 2.*F<  J» 

r ^ « 

'•  • 4 

73 

YU)  » SAYEYI  J)  a.S»H*FU) 

'“3  1 <- 

RUNGE«t 

^ A 

RF7U  -N 

'•T  ’ 1 

«• 

DC  HA  J«l.N 

^ • 4 -» 

FHI  ( J)  > 7PHI  ( J)  ♦?,  •FIJI 

AA 

vlj 1 ■ SAVFYI J)  ♦ HAF(J) 

«• 

Y ■ Y ♦ r. 

^ 4 - 

eu\‘.  = » 1 

R"“U?  s 

*"1  "• 

DC  3S  J-1,N 

■:  ^ 

•’(Jl  » SAVFYUlAl  P“i  U)aF(  1)  )•M/^. 

'A  -A  -; 

RUNG-  > « 

:-Nr 

* 

• ■ I\  rPr 

rr’-  n.£HC&rc. scukcf .noi is". v?’'FC'< .i -a:  - 

NAf'E  » » iiN=c\’‘  • 

*::s*  SCUOCE  STATf^fNTS  - 1?.P5r:,KA''  <^17!  * 

■•C!?*  lAiiNTS"  IC  ' '■.FNERATcrt 


IV 


Hf  A^rv:  co**®iLA''ioNS 


App«ndlx  B 


Equations  for  optimal  control  tasad  on  the  condition  that 
relative  suspension  mean  deflection  is  1/3  of  allowable 
axle  clearance  H and  absolute  acceleration  transmissibility 
be  a minimum} 

for  n^^  harmonic  (j“  |P?) 

mj^5lj^+c(Xjj-x^)+k(x^-x^)-0  a) 

m^x^-c(Xjj-x^)-k(Xjj-x^)-*-Kx^-Kh„  exp(ju)nt)  b) 


Xjj-Xa-y  c) 

Multiply  a)  by  m^,  b)  by  nij^  and  subtract  b)  from  a)  and 
divide  by  } 

y^c(  (m^+mj^)/ro^mj^)y+k(  (m^+mj^)y+y  {K/m^)-XjjK/m|^)- 

^ (K/mjj)exp(  juu)t)  d) 


and  from  a) 

<c/mj^)+y  (k/mjj)-  0 e) 


Let  ®xp  j(nwt+^),  y"Y  exp  j (nuit+ip) 

Then 

(-u)^n^+  k-t-K  + k + c ( (m^+mj^) /m^n^)  jwn)  f -K/m^ 


m 


j“n  + k ) 


2 2 

ui  n 


f) 
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exp(  j«Ji)Y 


5in 


i) 


exp(jij»)X  - 


Kh^{cja)n/m^  + 

6 ra_ 


j) 


After  Rationalizing  i),  j) 


n 


n 


k) 


“a  fia 


l+4c^n^ 
. ni. 


1) 


-2 

y 


t i. 

I\  % 

7 

a 

7 Z 

- On  (1+k  + % 

' K “? 

K 

A 

Z 2 

)+ll  + 4i:  On  In 

in  k)  , 1) 

2 

n * wn/wj^  , 

, K/m^ 

c/Cc  , 

- 2{k  ; 

: be  required  that 

< H ot 

y = H/a  - H/3 

(a-3) 

- 1 

? 21S 

m 2 

S h2  (X/h  )'' 

■*•  n 

= 1 

2if 

ITt  A 

z y^ 

1 

OK  ~T 
n 

A 


m) 

n) 


x' 

2-2  -A 

- (l+4c^y‘') 

-2 

n 

> mean  square 

n2 

■ /(H^  ” 4ttu^ 

= mean  square 

The 

program  for  n)  and 

o) 

follows. 

2 


speed  u 


P) 


62 


Mt 


MAIN 


MIC^Ii'AN  TfRujNAi  SYSTE"  FOtTRAW  C14niM 


0^-0l-7<> 


THE  VEHICLE  tf^AKAHETERS 

AM  - BOOf-HAl^-AM  « AXLE  MASS*  TK  - TIRE  SFPIN' 
RATE*  SX,;*  IttifENSION  SPRING  RATE*  C - DAMPING 
COEFF  icrtHtffiii  . number  of  data  POINTS*  V • SP 


IFOR  FEI 


lx*  WHERE  N IS  AN  INTEGER} 


0001 

0002 

0003 

OOOA 

0005 
0C06 
0007 

0006 

0009 

0010 
0011 
0012 

0013 

0014 

0015 
CO  16 

0017 

0018 

0019 

0020 
0021 
0022 


0023 


implicit  RCML  PSIA-H*!)-!) 

COMPLEX  * SHI  I 28  »•  SHNEWII28) 

DINENsfbN''i^TI2)*  WAItSI 

fyiMEAtSIfH  X#tt0l*Ati««^T}l8^i281v»X2M-12OV*0X2BI  126} 
01  MENS  ION  r(4l#DERt|4)*PRMTI5)*AUXI16»4) 
DINENSWN'1^128) 

oiMENsiBit  - 

PI  MENS  Ifli^JMfll»PHASEF<  1281 

coMMONyiWM||MHR  ^ 

COMPL  E X ^VO^Nl  1 29 ) • 0 1 F HI  1261 
DIMENSIOH  sm;2f}*  S21129)*  0PSDI1281 
C0MM0Ny¥(DIAm#AN*8M*TM*SK*CyPI  -- 

COMMON/XMBRltotf*  AL 
COMMON/RflNHlItfN- 

INTEGER  mvm 

PI  - 3.l4i9n939G9793 
R E A D I 5 » 24I|H1IIW8M » TK » S K *C 
REAOI5>mT»VfL>MCteAR*NP 

FORMAT  I mytaoT' 

FORMAT  I 

WRITEI6t«||fl[Mp»«l*^TK*SK*C~  . . .. 

WR1TEl6«Vi|[T*ilL*NP 

FORMATItlCl/^»«B«S  VEHICLE  SUSPENSION  DESIGN  %%%' i 


FORMAT  I mUBAt' 

FORMAT  I 

WR I T El  6«|Ptt|f»8M»TK *SK*C~ 
WR1TEl6«Vi|[T*ilL*NP 
FORMATItl||/W»*B«S  VEHICLE  S 
9 5X  * * 

35X*'  INPUT  DATA  FOR  THE  SYSTEM 
9 5X**  UNSPIlHlft  MASS 
;»5X**  SPWUWG  NASS 
35X.*  YitE  SPRING  Rate 
95X*'  SUSPENSION  SPRING  RATE 
a 5X  * * DA  MP  IMG  tOEFPlC  I ENT 

format  I9IC«>  liNPUT  0ATA:'< 


: '/ 

•*F10.2.*  KG'/ 
•*flO*2*^  KG'/ 
•*F10.2»*  N/**'/ 

**FI0.2**  N/M'r 
•*fl0.2»*  •/ 
'///  ) 


35X*' 

95X.' 


ISIIK*  IMPUT  data:'/ 
SP<W">R'  -the  vehicle 

TERfUlN  WAVE  length: 


.F7.  2, 


NUfWIR  OF  OBSERVATION  POINTS  - 


M^SEC  / 
*15// 

• I 


0025 

c 

U . 20.1^,  . 

.0026 

CALL  R0U^1NP»SH} 

0027 

00  12  l•l»12• 

0026 

12 

SSHI  1 1 -VSNI 1 1 

C 

00  500  |U»8»18 

0029 

00  30 

N • 7 

CALL  FFmil*>l*NPI 

\4-  , 

4 *>&%.•*  • * - . A- 


i \ 


• w'^**^**>*  V*  ** 


1 


O'"  ^ 1 

c 

CALL  0PTIH|AH»BH.TK»P1 »AL»U»SH,KP) 

003? 

CALL  PS0CAL«SH»S1*6P) 

OP33 

CALL  FUNIFI 

0034 

DO  14  — 

CC35 

14 

PHASEFU)  ■ Fin 

00  36 

00  19  I«2*l|t 

0037 

19 

Fin  ^mniiirn  — 

0030 

FMI  • -mattFIlZSI 

0030 

T -0.0 

004C 

t 

C 

TMAX  • AtFV  ' 

0041 

OELT  • - - - - 

004? 

C 

NOIH  - 4 

CC4  3 

c 

DC  ? I-l»N0|W 

0044 

? 

YUI  - 0.0 

0045 

3 

K • R w NOC  o w m * • Fi’  w V C»  T I • 

0046 

IFlK.NE.lliO  TO  9 

0C47 

OERril)  - fl3> 

0048 

cePY<?i  - vw-  - 

0040 

CCC  » C*iyi4»  - YI3n  ♦ SK*IY*2>  - Till) 

00  50 

X - TFDEIT 

0C51 

n - IF|X49flOtr4KH 

005? 

PESO  ■ X-II 

005‘« 

IF(RESD.GT.0.79III«II4| 

005i. 

FF  - F«  Ilf 

0C55 

XNEW  -II 

0056 

IFIDABSIX-XNEWl.LT.O.lICO  TO  7 

3057 

FF  • (PI  Il»|l»f  mil-FZ^O- 

0C5? 

GO  TO  6 

0050 

7 

XA ( 1 1 ) - nil 

0060 

xBim  • y-w-F 

0061 

0X?4(m  • !l./Ani4ICCC-TK*I  Y(  n-FF)  » 

006? 

CX?BiIII  - 1 1.0FBHIP(-CCC» 

006^ 

6 

D£or(3)-M.O#A#»*ICCC-TK*(Y(  ll-FFM 

00c*> 

c 

DERYU)  - (l.0/BH)4(-CCC) 

006  5 

c 

5 

CONTINUE 

0066 

IP  (T.LT.THAXICD  TO  3 

r C2NSTBAINT  EOUATICS  fOP  QPTl HI ?aT | ON 


C 

c 

oct? 

301 

FORHATI •ROOTS  ARE*#  2E17.7) 

00o8 

V 

CALL  HCAL-(0X2A»0A2B#XA#HH#NF) 

0069 

CO  18  I-l#120 

007C 

18 

ORlGHt  11  - REAIISSHC  11  1 

0071 

00  10  I-1»HF 

007? 

10 

SHNEWIIl  - OCHFIXIHHI I 1# 0.002) 

0073 

CALL  FFrfSHNEM»N*NP) 

I 

t 

1 


: t '.A* 

"cRNlNAL 

00  ? 5 

0076 

CC7T 

0078 

t 

0079 

201 

0080 

202 

0081 

0082 

206 

0083 

C08A 

C 

0085 

0086 

00  87 
0088 

0089 

0090 

0091 

-- 

0092 

16 

OOQ’ 

210 

C 

C‘C94 

0095 

O^'Qb 

99 

call  PSOCALf SHNEM.S?»NP) 

WPITE<6»20II 
DO  1 J»l.NP 

OPSOf JI-IOPIGHI J))*10n.0/CRIGH( J ) 

WR  ir  E(  frrMZHlItlGtlt  J)  #HH(  J t »DPSOI  J ) 

FORMATIlHl/Zf#*  COHPARISION  BETWEEN  ORIGINAL  i COMPUTED  fiT 
35X»*TERIIAIN  DATA*  # lOX* 'C  OHPUTEO  0ATA*>10X»*X  C I FFE  PENC  E ' /f 

a-T  ■ 

F0RNArf5X»3E22.5f ) 

WRITEI6*206I 

FORNAH  teORT^HMK  -OUTPUT*^  

35X* • FREQUENCY* *8X» •PSD* •tAX» *1061  FRO  I *»9X»* LOG ( PSD ) ' ) 

ONEGA  - 2.0PPIPY/AL 

TAKING  LOGRITHNS  OF  PSO  AND  FREQUENCY  VALUES 
DO  16  I«l#NP 

PSO»SH  I »^TA^Al-/<  2.0*P1 1 
FRO  « I*2.0PPI/Al 
PSDL«OLOGIO(rSO) 

FRQL  « DtOSlOtfRQT  

IFIll.EO.D.OR.I  1.GT.69)  )G0  TO  16 
WR ITEI 6>  2101 FRO*PSO>  FRQL  »PSDL 

CONTINUE  

FCRNATI 5X#6El7.t» 


99  CONTINUE 
5CQ5  STOP 

J^'Qb  END 

-coticns  in  effect*  IO#EBCOIC»SOURC£»NOL ist»nodeck»load»nonap 
♦OPTIONS  IN  effect*  NANE  • NAIN  » LIMECNT  • 57 
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'•'M 
r, ''  ■>  T 

GC04 

0009 


OC  C 6 
000  7 

0000 

0010 


0011 

1 7 

. 

' I fc 


c 

c Suni;PUTINE  C^TI"  S^.*P7f  hE«5E 

r 

SUBROUTINE  CPT  I H ( AN.  B M»T  K » P I » At  .U»  SH,N«>» 

INPlICIT  REAl  *8IA-M.C-/I 
C 

OINENSION  R00TI91.  VIO) 

CIHENSION  C0EI9)>  RRI9I*  RCI9)»  PCLI9) 

COMPLEX  *16  CDNSH.  SHfl?8) 

C 

C WRITE! 6*A1I)U 

All  FORMAT!*  U -*.017.6) 

C WRITE!6.40n!SH!  1 1.1-129.128) 

SIGNAH  - O.H - 

C WR  ITE!6.A02)NP 

AC?  FORMAT! * NP  -* .19) 

00  10  N-I.NP 

10  SICMAM  « SICMAH  ♦ CDABS!SH!N)) 

C 

C WR ITE! 6.A07)SIGMAH  « 

A07  FORMAT! *SIGMAH  -*.017,6) 

C 

rO  AO  IH  » 19.  AO.  3 
XSO*IN  . 1.0019 
H * IH*0.01 
ANP  « NP 

CC  « 9.0«SIGMAH/( ANP-M-H) 

WR  I'^E!  ^.«>oe)cc 

PCFMATI'CC  -*.017.6) 

CC  3C  17  - 1.10 
Z « 0.1*IZ 
WR 1TE(6»A21) 

A?l  FORMAT!*  ENTEREC  FIRST  OC  LGCP*I 
CC  30  ix.  1,10 
5K  - 1 000.0* IK 
CMB  « CSORTISKFBM) 

C*'A  , rSORT(TKZAM) 

CMR  * n*13FGMA 
WR  I-^E!  0.-22) 

>,->2  FCPMatc  entered  ■'HE  SECCND  LCCP*) 

SKR  . SK/TK 
n 1 * SKR  ♦ OMRPQMR 

C12  r 1 1 ♦ i.o 

»2  . 2-7 

r-Ed)  * CC 

CCE!2I  ^ C.O 

rofc!  3)  . ‘..•72  - ?.*C12 

COE(A)  • 0.0 

rCE!9)  - C12*C12  ♦ 2.*CMR*CMR  - CC  - 8 . • Z? -C * 1 -C  1 1 

cnE(6)  » 0.0 

rOE!7)  . A.*Z2*C  II*C  11  - 2.*0MR*0mr.cI? 

C0EI8)  « 0.0 
COE!'))  « CMR**A 
C 

c 


66 


^.c  . 


4r?  FQRMATCCCE 
C CALLING  A PROGRAM  FROM  SS® 


C 

CALL  OPRRM|COE»9>fiR>RC*POL»iRp lEP) 

0045 

c 

404 

PORNATI*RRe  RC*>017.6) 

c 

WRITEI6p409IIR 

0046 

405 

PORMATCMR 

0047 

IPIIR.EO.OIGO  TO  30 

0048 

I - 0 - — - 

C049 

J-l 

0050 

0091 

19 

IFIRCf JI.EO.O.OOIIGO  TO  20 

IP  f J.xO.  IRtOu~TO  'TT 

0052 

J • J4l 

CC53 

GO  TO  19 

0054 

20 

IPtRRUl.GE.O.OODCO  TO  22 

0055 

IPf J.EO. IRIGC  TO  25 

0056 

J « J ♦! 

CC57 

GO  TO  19 

0058 

22 

I « I ♦ 1 * 

0059 

ROOT!  I 1 • RRU  1 

0060 

J • J 4 1 

0C61 

GO  TO  19 

25 

continue 

? 

IFII.EO.OIGO  TO  30 

C 

WRITE(6p201  )(  UpRCOTU)  )p  J>U  n 

'•r 

201 

FORMAT!*  ROOT!', 11, *1  - *,017.6) 

DO  27  J • 1*1 

27 

VIJ)  * 0,5PR0CTI J)*AL*3M8/PI 

DO  29  Jl«l,l 

■'r6“ 

IPIOA0S(U-V!21il.GT.O.2)GO  TO  29 

CO  28  J » 1,  1 

T07G 

RCC2  ■ ROCTI J)*ROCT! JI 

vv  * 1 

RG04  • «Ca2*RQQ2 

fN/' 

»SC  « (1.44.4T2*R0Q2)7IR004*CC ) 

^ 7 

IFIHSO.GE.XSOMINIGC  TO  26 

V 1. 

XSOMIN  - XSO 

•,  r » t 

SKCPT  » SK 

7CPT  . 7 

T G7T 

28 

CONTINUE 

wC  71 

WR  ITE! 6,202) <( J,V<J) ), J-l* 1) 

:r  7= 

202 

FORMAT!*  V!*,ll,*)  -*, 017.6) 

or 

WR ITE! 6*  203)SK,{,H 

- . 

203 

FORMAT!*  SK  -*,017.6,*  7 -*,DI7.6,’  « 

:.c;e: 

29 

CONTINUE 

00  8 3 

?0 

CONTINUE 

CC04 

WRITE!6,418)SK0PT,2QPT 

Ov;65 

418 

FORMAT! *SK0PT-*.01T.6**  ZOPT • • ,D 17 . 6 ) 

006O 

40 

continue 

0067 

RETURN 

OCflfc 

END 

*OPTiONS  IN  EFfECT*  1 0»  E«€01<: » SOURCEvNa  ISTp  NOOECK  >L  OAO>  NOnaP 
•OPTIONS  IN  EPPECT*  NAME  - OPTIH  * LINECNT  • 
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.'V 


2<;4i5 


1 


I 


[ 

I 


I 


f 


«> 


1 “ • 
- 

01"'’ 
':  JO* 
0005 
''OOh 
0007 
coje 
r'5''o 
''010 
''01  1 
'*01? 
'01  1 
''01  t, 
c 


* • t • • 

TJ»*»liriNf  ^’JNUI 
71179) 

cr*<riea 

ca^^fjN/rFf/i.NF 

: c«««3H/vJ94f##*'Hi)9«r«»lli,c*  f t 

-• 

or  « t«UN» 

7(1)  • 0.0 

tkhP  • 7lOVT|4i»l — - 

?P  5 J*1*l(» 

Ml  • 1*J 

ii  • 9H7II# 

4*'I  - *104*(N1  ) 

r«  . I ANirAS»)*ij) 

'?  • 5i*(irm 


roi 
OJl  ^ 
fOl  ■= 


rr  . C3S(Tm» 

tMI)  • 47(1)  ♦ l(  il*;ML«(Cr.SS  I 


''•>1  r 

w V.  V. 
^ y • 

- T*  1 
> (. 
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Appendix  c 


The  algebraic  formulation  of  the  problem  is  given  next. 


rw  /r  K ^ 1 

(3) 

- 2 ,„n,2,l/2 

iMyl  ) 

(4) 

-2 

where  h = h .h*/2 
n n n 

(5) 

* th 

is  the  complex  conjugate  of  h^^  the  amplitude  of  the  n 

order  harmonic  component  of  the  terrain  elevation  spectrum  h^  . 


2 


(6) 


whore 

„ ^,2  , 2 2 2^.2, 

N^=K  (us  c n +k  ) (7) 

n,  I k (1- (n  ' )^(m  n^ui^-(k+K)  )+k^]^+  (8) 

i O fi 

2 *^2 
( ('n)“|i  i (M  ) -K+k  (ni.'/ii-  ) *■  1 

O 

|i-l”  N2  (9) 

where 


= 1 (10' 

* 

:■  -[(..>  .-n  ^ “^-1+ (n-.O  (!->.■  /n  . ) ^ ) /K- (k/K)  1 " + 

o a o 

I . cn)  (m^..^+k)  /K)  " 'k^  (11 

use  tile  fol  lowi  nq  nomenclature: 

: eircular  frequency  of  spruna  mass  when  stand i!i<' 

> •.  .int  ' fad  'ser 

.ni'iessed  circular  ’'requency  rad  see 
: unsurunq  mass  spring  rate;  N/m 

k:  sprunn  mass  spring  rate;  N/m 

c:  damping  constant  of  shock  absorber  (average) Nsec/n 

. axle  mass  per  wheel  set; kg 

ti : integer  1,2,3, order  of  vibration 


75 


(24) 


B ,«n  u>  (c/ui^) 

X2  O 

Ayl-((a)Q/na))^-l+(n«)^.  (l-(Wo/na))^)ma/K-(k/K)  1 

(n^u)^  +k)/K 

Byi-1/K 


(25) 

(26) 

(27) 

(28) 


By2"(‘^cn/)c)  (Ayj/k+l/K) 

Substitution  of  (14)  to  (28)  into  (12)  and  (13) , respectively, 
will  yield  two  simultaneous  algebraic  equations  in  the  unknown 
parameters  k and  c.  Note  that  because  of  automatic  height 
control  with  equal  a constant  the  pro-rated  mass  m^  does  not 
enter  into  the  equations  explicitly  but  is  contained  in 

For  the  terrain  identification  (subroutine)  we  write: 

-2 


, n 

XX 


n n 
n 


R*  u/u)^"4i.  .u 
n n c 

-th 


(2‘>) 


where  in  (29)  is  the  n order  power  spectral  density  value 

(m^/rad)  , in  the  wave  number  of  the  n^^  order  wave  (Rad/m)  , 
n 

is  the  n^^order  circular  fraquancy  (rad/sac)  and 


n 


^ 'y 


th 


(m  sec/rad)  is  the  n order  power  spectral  density 
t n n ^ 


-•aluc  relative  to  time. 
. n 


So 


NN 


U 


(30i 


.1 1 so 


iiui 


Au  (approximately) 


t ” 


such  that 


(32) 


h > Au/.*’  - A/n  ■ ♦-w’O 
n XX 

The  average  A-  value  A^^^  is  then 


(33) 
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Partial  differentiating  (1)  makes  P a minimum,  if 
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ac  7c 
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^ y2 
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Aj^^-=k  (1-  (nw/w^)  (m^n^u)^-  (k+K) ) +k^ 

^x2*  (naj/u)^) 

®vi " ^ ) (K+m_n^a)^)  (K+m_n^u)^) +2k  (nui/to  ) 


avg 


n«100 
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(34) 


Equations  (1)  through  (34)  are  processed  and  the  program  is 
shown  next. 
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Computer,  Suspension,  Simulation,  Terrain,  Model,  Ride,  Vehicles,  Control, 
Spring,  Damper,  Ride  Control,  Adaptive  Control,  Random  Vibrations,  Power 
Spectral  Density  of  Terrains,  Alrsprlngs,  Hydropneunatlc  Springs 


This  refX)rt  contains  the  software  package  for  an  adaptive  optimal  suspension 
control  system  relative  to  terrain  random  vibration  dlstxurbances.  The  promsed 
problem  solution  is  shosm  to  fall  Into  two  separate  program  categories:  a) 
recognition  of  the  terrain  and  parameter  selection,  by  means  of  an  oi>-board 
minicomputer  or  microprocessor,  b)  optimlaatlon  of  suspension  parameters  for 
arbitrary  terrain  configurations  obtained  from  terrain  statistics  and  executed 
on  a centrally  located  statiotutry  computer  facility.  - — > 
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I^The  interface  bettMen  the  stationary  conputer  facility  and  the  on-board  micro- 
processor is  accomplished  by  means  of  a data  bank  prepared  at  the  stationary 
facility  and  peraanently  stored  in  the  memory  of  the  on-board  microprocessor. 

The  stMpenslon  parasmters  are  set  by  a servo-control  unit  on  the  vehicle  which 
is  activated  by  the  sdcroprocessor.  The  servo-control  'jnit  regulates  the  supply] 
and  release  of  air  in  ttie  hydropneumatic  suspension  system,  thereby  increasing 
or  decreasing  the  spring  rate  according  to  the  optimal  requirements.  In  a 
similar  ownner  the  damper  orifice  sise  is  increased  or  diminished  depending  on 
the  required  effective  damping  parameter.  If  need  arises,  the  vehicle  can  opex^ 
ate  at  fixed  suspension  parameters.  The  results  of  the  investigation  are  shown 
in  graph  fora. 
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